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THE INFLUENCE OF VERTICAL WIND GRADIENTS ON THE 
LONGITUDINAL MOTION OF AIRPLANES 

By Joseph Gera 
Langley Research Center 

SUMMARY 

The present study is an attempt to make an assessment of the influence of wind 
shear on the longitudinal motion of airplanes. It was assumed that the wind is completely 
horizontal and its speed varies linearly with altitude. It is shown quantitatively that both 
glide and climb performance are influenced by wind shear and that trimmed flight at con- 
stant airspeed , attitude, and with fixed controls is along a parabolic path relative to the 
ground. The problem of the landing approach in a wind shear is examined in some detail. 
Small -disturbance theory indicates no wind-shear effect on the short-period motion and 
the time for the phugoid to damp to half amplitude but the phugoid frequency and damping 
ratio vary considerably with wind shear. A nondimensional quantity which depends on the 
wind shear and airspeed is shown to be a fundamental parameter influencing the longi- 
tudinal dynamic behavior of the airplane. 

INTRODUCTION 

All-weather operation on a routine basis in civil aviation is one of the major goals 
of the aviation industry. The most critical phase of all-weather operation is the landing, 
as attested to by the establishment of the various weather-minimum categories for landing. 
While progress is being made in the design of ground-based and airborne landing aids, 
adverse weather conditions — especially variations in wind direction, strength, and gusti- 
ness — remain a major factor affecting safety during landing (refs. 1 and 2). In the con- 
sideration of wind effects a distinction should be made between random gusts and a veloc- 
ity gradient in the mean wind or wind shear. Random-gust or turbulence effects are 
important with respect to the structural design and the riding qualities of the airplane. 

In contrast, wind-shear effects become significant when the airplane is to be controlled 
with a high degree of precision as, for instance, during an automatic landing. According 
to Henley (ref. 2) , wind shear was the main cause of touchdown error and of the large 
variations in rate of descent during the evaluation of a particular automatic landing sys- 
tem. It seems probable that other landing systems will be similarly affected unless 
control-system designers take wind-shear effects into consideration during the selection 
of the control laws for their system. 



Although the effects of wind shear on the dynamics of the airplane have been con- 
sidered in numerous engineering analyses, experimental data are scarce. Some early 
NACA flight tests reported in reference 3 conclude that landing approaches made by con- 
ventional airplanes while deliberately disregarding wind effects are unduly hazardous. 
Reference 3 also includes actual measurements of wind shear near the ground. Refer- 
ence 4 presents the frequency of wind-shear occurrence and magnitude in a particular 
locality along with the results of a flight test showing that the effect of wind shear on the 
climb performance of the B-29 aircraft was significant. 

Analytical studies of the effect of wind shear on the longitudinal dynamics include 
references 5 and 6. Additional information including the effects of other atmospheric 
gradients such as density and temperature variations can be found in references 7 and 8. 
In some relatively recent publications there is a lack of agreement as to what actually 
takes place during an approach to landing in a headwind which is decreasing with altitude. 
According to references 9 and 10 an airplane trimmed for a glide at constant airspeed and 
constant attitude would touch down short of the initial aim point if no control action is 
taken by the pilot. In contrast, reference 11 asserts that a trimmed airplane executing 
an instrument approach in decreasing headwind would overshoot the intended point of 
touchdown unless the pilot reduces power. 

The present study is a simple attempt to make an assessment of the influence of 
wind shear on the longitudinal dynamics of airplanes. Improved knowledge of wind-shear 
effects is important in the future design of approach and landing aids. 
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gravitational acceleration, 9.80665 m/sec^ 



i,j,k orthogonal unit vectors 

<v IW 
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rate of change of pitch angle, rad/sec 
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time, sec 
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velocity relative to air, m/sec 
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perturbation velocity, m/sec 
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u’ 
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horizontal wind velocity, m/sec 

, du,„ _i 

wind velocity gradient, , sec 


~ earth 
X,Z 


inertial velocity, tailwind positive, m/sec 

aerodynamic and propulsive force components in the i and 1 
respectively, N 
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linear displacements in earth-fixed axes, m 
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direction, 
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‘-'OL 
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rad-sec^ 
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1 

az 

m 
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m 

ad’ 

rad-sec 

z 6 = 

1 

az 
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m 

36’ 

rad-sec^ 


a angle of attack, rad 

/3 exponent in the wind-velocity relation 

r angle between earth-fixed x-axis and the unit vector j., rad 

y perturbation from the equilibrium value of r, rad 

y Q flight -path angle seen by an observer moving with steady wind, rad 

6 generalized control displacement, rad 

£ damping ratio 

0 pitch angle of the airplane relative to earth-fixed axes, rad 

8 pertubation from equilibrium value of 0, rad 

p atmospheric density, kg/m^ 

, Un u W 

a nondimensional wind- shear parameter, 

g 

co circular frequency, rad/sec 

Subscript: 

o reference condition 

A dot denotes differentiation with respect to time. 
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ANALYSIS 


Wind Shear 

In the technical sense the term wind shear means a spatial variation in the mean 
wind speed. According to reference 12, even in the presence of mechanically induced 
turbulence the mean-wind profile is adequately represented by the law 

0 


where (u w ) ^ and (u w )g are mean wind speeds at heights and z 2 , and 0 usually 
takes a value of 1/7 in an open terrain. In this report it is assumed that 0 equals 1 , 
so that the mean wind speed varies linearly with altitude according to the simple law 


u^ being the wind velocity gradient whose algebraic sign depends on the orientation of 
the airplane in the moving atmosphere. Figure 1 is a sketch which illustrates what is 
meant by positive and negative wind shear. Absolute values of u^ up to 0.03 sec - * 
are fairly common in the lower layers of the atmosphere with extreme values reaching 
0.15 sec - *. 


H (*J\ 

( u wj 2 \ z 2j 


Altitude 

(-z) 



Negative shear, u' w <0 


Altitude 

(-z) 



Wind velocity 
Positive shear, u' w >0 


Figure 1.- Definition of positive and negative wind shear* 
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Equations of Motion 


This investigation is concerned with the effects of wind shear on longitudinal dynam- 
ics of airplanes; consequently symmetrical flight of the airplane as a rigid body will be 
assumed. Figure 2 shows the coordinate systems used in the analysis. 




Figure 2.- 


Coordinate systems. 


It is assumed that the earth-fixed axes are fixed in space. The moving system of 
unit vectors { i.i.k') are defined as follows: the unit vector i points in the direction of 
the relative wind, the unit vector k is directed opposite of the lift vector, and i is 
defined by the requirement that (i, j,k) be a right-handed system. The absolute angular 
velocity of the moving coordinate system is seen to be T j. 

The absolute linear acceleration expressed in the moving coordinate system is 


Y ea rth = 3F( U + Uw C ° S F )~ + 5l( Uw Sin ^ + ^2 X [( U + u w cos r )i + u w sin 

where x denotes vector cross product. This expression can be simplified to 

dearth = ^ + i w cos ^ + (“ u r + “w sin r )iS 
The time derivative of u w may be eliminated from this expression by noting that 

i » = 4 u »ab = ‘V ,sinr 

so that the final form of the linear acceleration becomes 

V = (U - Uu’ sin T cos r)i + f-Uf - Uu’ sin 2 r)k 
-earth V w V w J~ 

Figure 3 illustrates the external forces acting on the airplane. All aerodynamic 
and propulsive forces are included in the X and Z components. As is shown, the 
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k 

/v 


I 

Z mg 

Figure 3«~ Forces acting on the airplane. 

origin of the moving coordinate system coincides with the center of mass of the airplane. 
The complete nonlinear longitudinal equations of motion are 

m(u - Uu^ sin T cos r + g sin r) = X 
m(-Ur - Uu^ sin 2 r - g cos r) = Z 
I y 0 = M 

where 0 is the angle of pitch and M is the total external moment acting on the air- 
plane. It should be noted that for zero wind shear the above equations become the stan- 
dard equations of longitudinal motion referenced to the wind-axis system. 

Equilibrium Flight Path 

As is usual in airplane stability and control analysis, it is assumed that the motion 
of the airplane consists of small perturbations from a reference condition. In the pres- 
ent report the reference condition denotes flight at zero external moment and constant 
airspeed, angle of attack, and pitch angle. It should be noted that according to this defi- 
nition of reference condition, the airplane may be accelerating along a curvilinear path 
when viewed by an inertial observer. Denoting the reference condition by the sub- 
script o, the equations of motion become 

m(-U 0 u^ sin T Q cos r o + g sin T 0 ) = 

m (" U o u w sin 2 r o - g cos Tq) = Z 0 \ (1) 

0 = Mo J 
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( 2 ) 


Thus, on the reference path one obtains the condition 


tan T n = - — — — 7 

Z Q + mU G u’ 


w 


Equation (2) may be used to show the influence of wind shear on glide or climb perfor- 
mance on the reference flight path by relating r Q to the flight-path angle y 0 that the 
airplane would have in still air at the same airspeed, angle of attack, and power setting. 
With the assumption that the atmospheric shear has no influence on the local flow around 
the various components of the airplane, equation (2) can be rewritten as 

_j \ tan ^o 


tan r n = - • 


X r 


il + 


intJqUw 


1 - a sec 


7 


^o 


(3) 


where a is the nondimensional quantity \J 0 \i'/g. In the remainder of this report 


will be referred to as the wind-shear parameter. 

<r 


Figure 4 is a plot of sin r Q against 


t? 


CO 



Sin y« 


Figure 4 '.- Effect of wind shear on the nondimensional rate of climb and sink. 
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sin y Q for various values of a. The sine of the respective angles is chosen since it is 
a direct measure of the normalized rate of climb or sink: 

y Q (no wind) 

T 0 (wind shear) 

Equation (3) shows that for positive values of a the angle T 0 is greater in absolute 
value than y . Thus, an airplane which climbs into an increasing headwind or decreasing 
tailwind experiences a rate of climb higher than the value it could achieve in perfectly 
still or uniformly moving air. Negative values of a would be reflected in an opposite 
effect. This phenomenon is well known and the utilization of wind gradients by certain 
kinds of birds is called dynamic soaring. Wind-shear effects in airplane climb perfor- 
mance were observed and measured during the flight tests reported in reference 4. From 
the wind-velocity measurements which accompanied the flight tests of reference 4, the 
value of a may be calculated. Based on the true airspeed of the test airplane the value 
of a is calculated to be ±0.105, the algebraic sign depending on the direction of flight 
relative to the wind. With the known value of a, the net increase or decrease of climb 
rate due to wind shear may be calculated from equation (3). The results of this calcu- 
lation are presented in figure 5 along with two flight -test data points taken from refer- 
ence 4. Instead of using the exact calculations in equation (3), reference 4 proposed that 
the net increase or decrease of climb rate due to shear may be obtained approximately 
by multiplying the rate of climb in still air with the value of a. For small climb rates 
this is a reasonable approximation as shown by the broken lines in figure 5; for large 
rates of climb, however, the use of equation (3) is recommended. 

Examining the equilibrium flight path in wind shear, one finds that it is not recti- 
linear in the earth-fixed coordinate system. This can be seen by integrating and com- 
paring the equations 


/- Uo sin y Q 

(steady wind) 

\_-U 0 sin T 0 

(wind shear) 

fu o cos y D - u w 

(steady wind) 

\u 0 cos r 0 - U w (z) 

(wind shear) 


in which a steady headwind of magnitude u w is assumed and in the case of wind shear, 
the dependence of the headwind on altitude is expressed by writing u w (z). It is inter- 
esting to compare the trajectory of an airplane in a steady glide into a uniform headwind 
with the flight path of the same airplane flown in a wind shear. In order to make the 
comparison more meaningful, let it be supposed that the initial height above the ground 
is Az and that at this height there is a headwind of velocity u^ Az in both cases, 
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but in the case of wind shear the wind diminishes to a no -wind condition at ground 
level. With these assumptions the equation for horizontal velocity becomes 


U Q cos Yp ~ 

U 0 cos r Q + u^(z - z Q ) - u^ Az 


(steady wind) 
(wind shear) 


This pair of equations can readily be integrated along with the equations for z, so 
the horizontal displacement x can be evaluated at both conditions for the same height 
loss. After some manipulation the following expression is obtained: 


Ax 


- «shear ’ «*eady wind = W 


Making use of equation (3) results in the simpler form 


Ax = 


sin 


Vo 


Az + 


g 


/sin y 


2 (Az)' 


2U 0 2 \ sin r o 

This expression gives the amount of overshoot or undershoot at ground level relative to 
the flight path in steady wind. Thus, in a headwind decreasing with altitude there is an 
undershoot (Ax < 0) as long as the inequality 

2U~ 2 

Az < 


g 2 - 


sin y 0 
sin r o 


is true. If the height loss during the descent through the shear layer exceeds the right- 
hand side of the inequality, an overshoot will occur. It is important to note, however, 
that for moderate values of wind shear there is always an undershoot initially regard- 
less of the thickness of the shear layer. 

In the above analysis the undershoot was calculated for the case in which the air- 
speed, angle of attack, and throttle setting were the same in both the steady wind and the 
wind shear. It was found that these conditions were possible if the airplane in wind shear 
assumed a pitch attitude different from the steady-wind value. 

Another comparison is possible by requiring that the airspeed, angle of attack, and 
pitch angle be the same for both the conditions of steady wind and wind shear. This is 
possible, at least in theory, by having different control deflections and throttle settings 
under the two conditions. Since now T 0 = y Q , equation (4) can be simplified to give 

Ax = lx) - (x) - S a (Az ) 2 

w shear v ' steady wind 2 sin y 

Li Uq O 

This value of Ax in a normal glide is always positive, hence an overshoot will occur. 
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Height, meters 


Figure 6 is a summary of these results for the case of a hypothetical STOL airplane. 
Thus, the question whether an overshoot or undershoot occurs can be answered by ana- 
lyzing the method the pilot chooses in making an approach in wind shear. Relative to the 
steady -wind case, an undershoot will occur if in a wind shear the airspeed, angle of attack, 
and throttle setting are kept the same. If on the other hand, the pilot changes the throttle 
and control settings so that the airspeed, angle of attack, and pitch attitude are the same 
as in the steady wind, an overshoot takes place. 

Qualitatively these conclusions remain valid for wind shears with other than the zero 
surface wind considered here. 



Figure 6.- Approach profiles of a STOL airplane. Controls fixed; airspeed = 36.04 m/sec; 

wind shear = 0.0845 sec-1. 
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Wind-Shear Effects on Longitudinal Stability 

The assumption of constant airspeed, constant angle of attack, and pitch angle is 
basic to the concept of trimmed flight. In the preceding analysis it was found that these 
conditions of trimmed flight in wind shear resulted in a flight path which was not recti- 
linear for nonzero r Q but appeared as a parabolic segment when viewed by an observer 
on the ground. (It should be pointed out that for the pilot without accurate attitude refer- 
ence the curvilinear path would appear as the normal condition of equilibrium flight.) The 
tendency of an airplane to remain on the curvilinear flight path in the presence of distur- 
bances is now examined by the aid of the usual small perturbation theory of airplane sta- 
bility. Thus, it is assumed that U = U 0 + u, r = r 0 + y, and © = © 0 + 6 in the equa- 
tions of motion and that the perturbation angles y and 6 are small along with the 
products of all perturbation quantities. If the aerodynamic and propulsive forces and 
moments are expanded in the usual manner and the substitution a = 9 - y is employed 
the linearized equations of motion are obtained as follows: 


d 

dt 

-|u^ sin 2r o - X u 

-X a 

g(cos r o - a cos 2r 0 ) 

u 


X 5 


-Zu - ^ sin 2 r 0 

"( z a + ' z “ 

-( u o + z q)^ + g sin r ° ' u ° u w sin 2r o 

a 

= 

Z 6 


-M u 


s(s - Mq) 

y 


m 6 


Stability derivatives appearing in the matrix equation are standard and their defini- 
tions are given in the list of symbols. The derivation of the above equation is straight- 
forward, and since it is done entirely within the framework of conventional airplane sta- 
bility theory it is not given here. 

On the right-hand side all control-dependent terms are grouped together. There is 
seldom a need for the complete form of the above equations; usually the rotary deriva- 
tives and Zq are small enough to be neglected. The perturbed motion can there- 

fore be characterized by the stability determinant 

s -|uV sin 2r 0 - X U -X a g(cos T Q - a cos 2T 0 ) 

-Z u - u^ sin 2 r 0 -Zq, -U Q s + g(sin r Q - a sin 2r o ) 

-M u s 2 - + Mq^s - Mq, s(s - Mq) 
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where s is the complex variable. In case of no wind shear the determinant reduces to 
the conventional longitudinal stability determinant written in terms of the motion variables 
u, at, and y. 

It will be assumed that even in the presence of wind shear the perturbed longitudi- 
nal motion consists of the usual normal modes; i.e., a lightly damped phugoid mode and a 
well -damped short -period mode. In the latter mode airspeed changes are considered 
negligible and the motion is approximately characterized by the determinant 

Z a U Q s - g(sin r G - o sin 2r 0 ) 

= 0 (5) 

s2 ‘ ( M <i + M< l) S " M « s (s - Mq) 

Thus, the presence of wind shear is reflected only in the gravity term in the upper 
right-hand corner. Since gravity effects are not of great importance in considering the 
short -period mode regardless of the orientation of the gravity vector, it can be concluded 
that the influence of wind shear on the short -period mode is a minor one in comparison to 
the influence of the derivatives M^, Mq, and Mq,. 

The long -period or phugoid mode is approximated by assuming that the inertial and 
aerodynamic damping contributions to the total applied moment are much smaller than 
those due to changes in u and a. The approximate stability determinant of the phugoid 
mode is 

s - | u^ sin 2T 0 - X u -Xq, g(cos T 0 - a cos 2r 0 ) 

-Z u - u^ sin 2 r 0 -Zq, -U 0 s + g(sin T 0 - a sin 2T 0 ) = 0 

— — Mq; 0 

Or, in an expanded form the determinant may be written as 

s 2 + 2£ p o> p s + uj p 2 = 0 (6) 

where 


Xq, g sin r o 
2 ?P“P = - x - + M u if' 
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The time required for the phugoid to damp to half amplitude increases in a climb 
and decreases in a dive as expected. Positive wind shear (a > 0) is seen to amplify this 

p 

effect. For sufficiently large positive values of the wind-shear parameter, Wp* becomes 
negative, regardless of the, sign of r 0 , and the subsequent motion is characterized by two 
real roots of equation (6). One of these roots can become positive indicating a diverging 
aperiodic mode. It is interesting to note that wind shear affects the phugoid mode even 
in level flight (r 0 = 0). The level flight case was used to verify the assumption of the two 
normal modes for atypical light aircraft whose aerodynamic characteristics and param- 
eters appear below: 


X u = -0.0451 sec 


-1 


Z u = -0.3697 sec - * 
M u = 0 5T3ec 


X a = -7.8729 
Z a = -108.79 


m 


rad-sec^ 

m 

p” 

rad-sec' 1 


Mq, = -8.811 i 

rad-sec^ 

M • = -0.90904 1 

a rad-sec 


“d = -2.0767 jjji 


sec 


Uq = 53.64 m/sec 


T 0 = 0 rad 

This set of data was used to solve for the roots of the characteristic equation asso- 
ciated with the above 3x3 determinant as well as for the roots of the approximate equa- 
tions (eq. (5)) for the short period and phugoid modes while allowing a to vary between 
-2 and +2. The short -period roots were found to be approximately -2.51 ± 2.60i with 
negligible variations due to changes in a. The short period approximation yielded roots 
very close to the exact values. The location of the exact and approximate phugoid roots 
are shown in figure 7. For values of a greater than 1 an unstable root appears; for the 
particular airplane this would happen if the wind shear exceeded the rather unlikely value 
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Figure 7.- Root locus plot showing the effect of wind shear on 

the phugoid roots. 

of 0.183 sec - *. It should be noted, however, that for higher reference velocities or at 
positive values of r o phugoid roots would appear at lower values of For certain 

configurations and flight conditions the appearance of this diverging first-order root may 
be important and warrant an analysis. 
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CONCLUDING REMARKS 


A study made to assess the influence of wind shear on the longitudinal dynamics of 
airplanes found that a linear variation of wind with altitude affects the longitudinal motion 
of the airplane in two ways. First, the path flown at constant airspeed, pitch angle, and 
fixed controls is not rectilinear with respect to the ground but rather it is a parabolic 
path. Second, the longitudinal stability characteristics of the airplane about the refer- 
ence conditions vary as a function of wind shear. 

Associated with the first effect is the apparent ability of the airplane to climb faster 
in an increasing headwind in comparison to flight in still air at the same airspeed, angle 
of attack, and power setting. Similarly, glide in decreasing headwind results in a higher 
rate of sink than what would be experienced in uniform wind or still air. If an airplane 
is trimmed during an approach in headwind decreasing linearly with altitude, it would 
normally touch down short of the point it would reach if the wind had remained at its ini- 
tial value. It is shown in the report that an overshoot could occur if the shear layer were 
thick enough. An overshoot takes place, also, regardless of the thickness of the shear 
layer, if the pilot uses the same airspeed, angle of attack, and pitch attitude that he would 
normally establish for a steady- or no-wind condition. It follows from these results that 
some elevator and throttle activity would be necessary to "rectify” the curvilinear path of 
trimmed flight as would be required, for instance, on the glide slope of an instrument 
landing system. 

Longitudinal stability of the airplane is influenced by wind shear even in level flight. 
The results of the investigation show that while wind shear has negligible effects on the 
short -period mode and the phugoid damping, the frequency of the phugoid motion and the 
damping ratio vary considerably with wind shear. In the situation of increasing headwind 
or decreasing tailwind with altitude the phugoid mode may become unstable if the non- 
dimensional wind-shear parameter is sufficiently large. 

Langley Research Center, 

National Aeronautics and Space Administration, 

Hampton, Va., August 9, 1971. 
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